Exact simulation of the Wright-Fisher diffusion
نویسندگان
چکیده
The Wright-Fisher family of diffusion processes is a class of evolutionary models widely used in population genetics, with applications also in finance and Bayesian statistics. Simulation and inference from these diffusions is therefore of widespread interest. However, simulating a Wright-Fisher diffusion is difficult because there is no known closed-form formula for its transition function. In this article we demonstrate that it is in fact possible to simulate exactly from the scalar Wright-Fisher diffusion with general drift, extending ideas based on retrospective simulation. Our key idea is to exploit an eigenfunction expansion representation of the transition function. This approach also yields methods for exact simulation from several processes related to the Wright-Fisher diffusion: (i) its moment dual, the ancestral process of an infinite-leaf Kingman coalescent tree; (ii) its infinite-dimensional counterpart, the Fleming-Viot process; and (iii) its bridges. Finally, we illustrate our method with an application to an evolutionary model for mutation and diploid selection. We believe our new perspective on diffusion simulation holds promise for other models admitting a transition eigenfunction expansion.
منابع مشابه
Markov Chain Monte Carlo Simulation of the Wright-Fisher Diffusion
Markov Chain Monte Carlo Simulation of the Wright-Fisher Diffusion
متن کاملWright–Fisher exact solver (WFES): scalable analysis of population genetic models without simulation or diffusion theory
Motivation The simplifying assumptions that are used widely in theoretical population genetics may not always be appropriate for empirical population genetics. General computational approaches that do not require the assumptions of classical theory are therefore quite desirable. One such general approach is provided by the theory of absorbing Markov chains, which can be used to obtain exact res...
متن کاملEfficient Forward Simulation of Fisher-Wright Populations with Stochastic Population Size and Neutral Single Step Mutations
In both population genetics and forensic genetics it is important to know how haplotypes are distributed in a population. Simulation of population dynamics helps facilitating research on the distribution of haplotypes. In forensic genetics, the haplotypes can for example consist of lineage markers such as short tandem repeat loci on the Y chromosome (Y-STR). A dominating model for describing po...
متن کاملRenormalization analysis of catalytic Wright-Fisher diffusions
Recently, several authors have studied maps where a function, describing the local diffusion matrix of a diffusion process with a linear drift towards an attraction point, is mapped into the average of that function with respect to the unique invariant measure of the diffusion process, as a function of the attraction point. Such mappings arise in the analysis of infinite systems of diffusions i...
متن کاملOn the two oldest families for the Wright-Fisher process
We extend some of the results of Pfaffelhuber and Wakolbinger on the process of the most recent common ancestors in evolving coalescent by taking into account the size of one of the two oldest families or the oldest family which contains the immortal line of descent. For example we give an explicit formula for the Laplace transform of the extinction time for the Wright-Fisher diffusion. We give...
متن کامل